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Abstract 

It is shown that analyses at the electroweak scale can be significantly affected due 
to Landau pole effects in certain regions of the parameter space. This phenomenon arises 
due to a large magnification of errors of the input parameters mt, ag which have currently 
a 10 percent uncertainty in their determination. The influence of the Landau pole on 
the constraint that the scalar SUSY spectrum be free of tachyons is also investigated.lt is 
found that this constraint is very strong and eliminates a large portion of the parameter 
space. Under the above constraint the trilinear soft SUSY breaking term at the electroweak 
scale is found to lie in a restricted domain. 



1. Introduction 



Currently among the models beyond the standard model, only those with supergrav- 
ity [1] can accommodate the unification of the electroweak and strong coupling constants 
that is consistent with the LEP data [2], provide a natural mechanism to break super- 
symmetry via a hidden sector [3] and can achieve a radiative breaking of the electroweak 
symmetry breaking via renormalization group effects [4]. The minimal supergravity model 
(MSGM) is one of the most studied of such models [4]. The parameter space of the MSGM 
with the radiative symmetry breaking consists of only four parameters: 

mo, nig, A t , and tan/? (1-1) 

in addition to the top quark mass, m t (which is indicated by CDF to be m t = 174 ± 15 
GeV [5]) and the GUT parameters, Mq (the grand unification scale), and olq (the gauge 
coupling constant at Mq). Here in equation (1), mo is the universal scalar mass at Mg, m g 
and A t are the gluino mass and the trilinear scalar coupling at the electroweak breaking 
scale, M ew (we choose it to be the Z-boson mass, Mz = 91.187 GeV), and tan/? = 
(H 2 )/(Hi) is the Higgs vacuum expectation value (VEV) ratio at M ew , where (H 2 ) and 
(Hi) give masses to the up and down type quarks, respectively. The parameter space given 
by eq (1) has two branches, corresponding to the two signs of the supersymmetric Higgs 
mixing parameter, fi, since radiative symmetry breaking determines only the magnitude 
of this parameter. 

Although the MSGM has only four free parameters, in addition to the 19 for the 
Standard Model[and 137 + 19 for the minimal supersymmetric standard model (MSSM)], 
the parameter space is still quite large and one may put additional constraints on the 
model to further reduce the allowed range of parameters. Constraints that have been used 
in the past consist of proton stability for the SU(5) type model [6], the neutralino relic 
density constraint consistent with the current estimates of dark matter in the universe [7] , 
the b — > S7 decay constraint [8,9], as well as constraints on the lower bounds on the 
SUSY particle masses from CDF, DO and LEP experiments. Thus for models with proton 
decay, such as the SU(5) type models, proton stability requires mo > m g and tan/?<10. 
The assumption of R parity invariance implies that the lightest neutralino (Zi) is the 
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lightest supersymmetric particle (LSP) for almost the entire parameter space. It does 
not decay due to R parity conservation, and hence may be regarded as the candidate 
for the cold dark matter in the universe. Recent analyses have shown that there exists 
a large domain in the parameter space which satisfies both the proton stability and the 
relic density bounds [7]. One may also impose additional constraints on the theory such 
as the experimental constraint of the b — > s*y branching ratio. This decay is interesting as 
it is one of the few low-energy processes that is sensitive to physics beyond the standard 
model [10], and the first experimental measurement of this process has been made [11]. 
Recent analyses [8,9] show that there is a significant region of the parameter space where 
the branching ratio B(b — > S7) lies within the current experimental bound even when 
combined with the relic density and proton stability constraints. 

In addition to the constraints discussed in the previous paragraph, the internal con- 
sistency of the theory also puts stringent constraints on the parameter space. For example, 
the color SU(3)c group should remain unbroken [12] including the one loop corrections 
to the effective potential [13] when the radiative symmetry breaking is imposed, all scalar 
particles must be non-tachyonic, and the allowed parameter space should be such that 
theory remains in the perturbative domain.lt turns out that the restriction that the mass 
spectrum be free of tachyons imposes strong constraints on the parameter space of the 
model. An interesting aspect of minimal supergravity unification is its predictivity. There 
are 32 supersymmetric particles in the model whose massses are determined in terms of 
the 4 parameters of eq. (1.1). Thus the theory makes 28 predictions. However, there is 
an important issue regarding how accurately predictions of the mass spectra can be made. 
We shall show that the precision with which predictions are made can be significantly 
influenced by the proximity of the top quark to the fixed point in the top quark Yukawa 
coupling. It turns out that in certain parts of the parameter space, the proximity to the 
Landau pole can magnify errors in input data such as m t and «g (which are known to no 
more than 10 percent accuracy) by an order of magnitude. Thus in this region of the pa- 
rameter space an accurate prediction of certain subset of the SUSY mass spectra requires 
a knowledge of the input data to an accuracy which is an order of magnitude more severe 
than what is currently known. Further it is found that the proximity to the Landau pole 
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often leads to the lightest stop mass turning tachyonic which eliminates significant regions 
of the parameter space. 

The outline of this paper is as follows: In Sec. 2, we give a brief review of the Landau 
fixed point analysis for the top quark Yukawa coupling. We then discuss two cases in 
which one can approach the Landau pole. These cases correspond to fixing the input 
data (such as the value of the soft SUSY breaking trilinear coupling) either at the GUT 
scale or at the electroweak scale. In Sec. 3, we discuss the magnification of errors that can 
occur when one is close to the Landau fixed point. The effect of these magnifications on 
precision calculations at the electroweak scale is discussed. We also discuss in this section 
the breakdown of perturbation theory as we approach the Landau pole due appearance of 
log(e)/e terms (where e measures the nearness to the Landau pole). In Sec. 4, we discuss 
the condition that the scalar sector of the SUSY mass spectra not develop tachyonic modes. 
We then use this contraint to determine the allowed range of A t . Conclusions are given in 
Sec. 5. 

2. The Top Quark Landau Pole 

Since the radiative breaking of the electroweak symmetry involves a large extrapo- 
lation, from the grand unification scale Mq to the electroweak scale M ew: the stability 
of solutions to small corrections needs to be addressed. A number of possible corrections 
from various sources have been discussed in the literature. These include two and higher 
loop corrections to beta functions, threshold corrections at the GUT scale due to heavy 
particles resulting from the breaking of the GUT group, threshold corrections at the elec- 
troweak scale due to the light SUSY spectrum, loop corrections to the effective potential 
and corrections due to possible higher dimensional operators from quantum gravity effects. 
Some of these issues have been studied in [14,15]. All of these corrections have influence 
on precision calculations at the electroweak scale. In addition, there exist another type 
of corrections that may be magnified due to proximity to the quasi fixed point of the top 
quark Yukawa Coupling [16-20]. It is well known that the one-loop RG equation analysis 
using the gauge and Yukawa coupling evolution below Mq gives rise to a quasi-infrared 
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fixed point for the top quark Yukawa coupling corresponding to D — where [16]: 



Yq = Y ! , (2.1) 



with 



F(t) 

Bo = l-6Y t JJ (2.2) 

Here in equations (2.1) and (2.2), t = 2\n(M G /Q), Y t = h\ /(47r) 2 , h t is the top quark 
Yukawa coupling constant at the electroweak scale, Y is the value of Y t at the GUT scale 
(i.e. at t = 0) and Q is the running mass scale. E(t) and F(t) are two form factors defined 
in [16] which only depend on the gauge coupling RG equations: 

F(t) = [ E(t)dt, (2.3) 
Jo 

and at the one loop level, E(t) maybe solved as, 

E(t) = (1 + /? 3 )^(1 + + (2.4) 

where (3i = ai(0)6i/47r, 6j are the coefficients of the one loop beta functions, given by 
(&i, b 2 , 63) = (33/5, 1, —3), «i(0) are the gauge coupling constant at the GUT scale, 
CKj(O) = ac and ot\ = (5/ 3) ay with Y being the hypercharge in the Standard Model. 
From equation (2.2), one observes that a fixed point in the top quark Yukawa coupling 
corresponds to Y/ (t) = E(t)/6F(t). The existence of this fixed point implies a fixed point 
mass for the top quark: 

m{ = (8Tr/a 2 (t)) 1/2 (Y t f (t)) 1/2 M z cos9 w smf3, (2.5) 

where sin 2 9 W = 0.2328 is the weak mixing angle. Thus, for Mq = 10 16187 GeV, olq = 
1/24.0, and M ew = M z , one has E(t) = 12.906, F(t) = 263.954, and hence Y t f (t) = 
8.149 x 10" 3 and 

m{ = 2.163M z sin/3 = 197.25 sin/3. (2.6) 

Because of the coupled nature of the RG equations, the same fixed point appears in some 
of the other parameters of the theory. Near the fixed point small variations in some of the 



input quantities can be magnified in the output. Specifically, we shall show in the next 
section that the computed quantities in certain regions of the parameter space can be very- 
sensitive to variations in ao, mt and tan/? near a fixed point. 

As is evident from the discussion above all these sensitivities in the parameters can 
be traced back to the fact that in the parameter space one is very close to the top quark 
Landau pole, which is determined by D being close to zero. We note that D depends 
on the gauge coupling constant RG equations only through a combination E(t)/6F(t) and 
on Y t at the electroweak scale. It is thus important to investigate the influence on this 
quantity from variations of the input parameters such as the grand unification scale and 
the GUT coupling constants etc. We consider first the size of two loop effects on the 
computation of the SUSY mass spectrum. In the following we discuss the relative size of 
the two loop effects on ct's relative to variations in the value of ao- The two loop RG 
equations for the gauge couplings are 



where aii(t) = ctij(f)/47r, for i, j = 1, 2, 3. In our analysis we have used the following 
values for ai at the electroweak scale, 

«i = 0.016985 ± 0.000020, a 2 = 0.03358 ± 0.000011, a 3 = 0.118 ± 0.007. (2.8) 

Table 1 displays a fit to the three gauge coupling constants using Mg = io 16187 GeV and 
various values for ao- One concludes that ao = 24.1 gives the best fit to the numerical 
values for the three gauge coupling constants given in equation (2.8), both at the one loop 
level and at the two loop level. One now infers from Table 1 that the variation in the 
GUT coupling has much larger effects on the unification than the two loop RG effects! We 
can similarly analyse the two loop RG effects on the mass spectrum and specifically on 
the top quark Landau pole. Table 2 shows a comparison between the one loop and two 
loop evaluations for the form factors E(t), F(t) and their ratio. It is seen in Table 2 that, 
although the one loop and two loop calculations for E(t) and F(t) differ by several percent, 
their ratio, and thus the fixed point for the top quark Yukawa coupling constant, Y t J , differ 
only by less than a percent. Again, we observe that the changes in these quantities due 




da 



(2.7) 
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to changes in aa are much larger than the changes between the one loop and two loop 
calculations. 

Next we discuss the question of what quantities may become singular as we approach 
the Landau pole. This question depends critically on what is chosen as the input and 
what as the output. For illustrative purposes we shall consider two cases. In each of these 
cases we shall use mo, m 1( / 2 , tan/?, ag, Mg, and m t as input. However, in Case I we 
shall assume that the trilinear coupling A is specified at the GUT scale, while in Case II 
we shall assume that the trilinear coupling A t is specified at the elecrtroweak scale. We 
discuss each of these cases in some detail. We consider Case I first. Here the value of the 
trilinear coupling at the electroweak scale, i.e., A t is given by 



where the form factors, H 2 , H 3 are defined in [16] and below. From the above we see that 
for fixed Aq there is no pole in A t , rather one finds that A t is smooth as one approaches 
the Landau pole, and in this limit we find that A t takes on the value 



We see now that A t p is independent of the value of Aq and is determined in terms of other 
parameters. Thus there is a reduction by one in the number of parametrers in this case. 
Further, here none of the computed quantities, such as squark masses, etc., exhibit a pole. 
This case has been extensively discussed in [18]. For Case II, A t is assumed fixed, and A 
can be computed in terms of A t via the relation, A = Ar/D , so that 



A t = A Q D + (m 1/2 /m )(H 2 - 6Y t H 3 /E) 



(2.9) 



A tP = (m 1/2 /m ){H 2 - 6Y t H 3 /E) 



(2.10) 



A R = A t - (m 1/2 /m )(H 2 - 6H 3 Y t /E) 



(2.11a) 



where 



H 2 = (a G /47r)(— h 3 + 3h 2 + —h^ 
o lb 



(2.116) 



and 



h i = t(l + l3 i t) 



-1 




(2.11c) 



Thus A has a pole similar to the pole in the top Yukawa coupling constant. 
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We examine next the behavior of the Higgs parameter fx in the limit as one approaches 
the Landau pole. We begin by noting that in the radiative breaking of the electroweak 
symmetry, fj, is determined by fixing the Z boson mass Mz via the equation 



tan 2 - 1 2 



H 2 = H \ _ 2 f 2 - -Mf (2.12) 

where 



= TO o + (0-3/i + 1.5f 2 )a G mj /2 (2.13a) 
m# 2 = A mi/ 2 / - (fc^o _ ^ m o) + em i/2 (2.136) 

Here /i, /2, /, fc, a, and e are defined in [16]. To investigate the Landau pole limit we 
write the functions h, k, f, e in the form shown below using the identity Dq D = 1 where 
D = 1 + 6 Y t F(t). We have 

h = {3D - l)/2 (2.14a) 
k = 3Y t FD /E (2.146) 
f = -6H 3 Y t D /E (2.14c) 
e = (3/2)[(GiA, + Y t G 2 /E) + (H 2 D + 6H 4 Y t /E) 2 /3 + H 8 ] (2.14a 1 ) 
From the above we find that m 2 Hi is smooth while rrijj develops a pole, i.e., one has 

m 2 H2 (pole) = -W t (F/E)A 2 R /D (2.15) 

From eq. (2.12) one can obtain the limit of // 2 in the vicinity of the Landau pole to be 

V? = /j 2 r /D +/j 2 (NP) (2.16a) 

where \j? R is given by 

and where (i 2 (NP) stands for the nonpole piece of fx 2 . We see from the above that the 
limit of quantities depends strongly on whether A or A t are assumed fixed as we approach 
the fixed point. For the case when A is kept fixed the limit is smooth, while the case 



when A is kept fixed induces a pole in several other parameters in the theory. The effect 
of these poles on electroweak physics is investigated in the next section. 



3. Effect of the Landau Pole on Physics at the Electroweak Scale 

As discussed in Sec. II for the case where A t is fixed, several of the soft SUSY 
breaking parameters and the parameter fx 2 develop Landau poles. Computed quantities 
that depend on these parameters can also develop poles and show rapid variations as we 
approach the fixed point. There are several elements of the SUSY mass spectra which 
show such a behavior. We recall that the stop quark mass matrix takes the following form, 



rri~ L m t (A t + /icot/3) 

m t (A t + /jcot/3) m 2 



" ] (3.1) 

tR 



where 



2 1 1 8 12 

m h = 3 m o + 3 m k + ("T^ 1 + h + 3/3)^/2 + m 2 t + (- - - sin 2 6 W )M 2 cos 2(3 

(3.2) 

12 18 2 

m i R = 3 m o + 3 m k + (3/1 ~h + ^h)^Gm\ /2 + m\ + - sin 2 9 W M 2 cos 2/3. (3.3) 

and m 2 Hi and m 2 H are as defined in eqs. (2.13). In the following analysis we shall assume 
that A t at the electroweak scale is given and its value at the grand unification scale A is 
determined by eq. (2.11). 

The light stop mass is given by 



ra- 
ti 



1 r 



(m~ L +mf R ) - J(m 2 L - m 2 R ) 2 + 4( A t + fj, cot /3) 2 m 2 . (3.4) 



while m? 2 is given by a formula identical to eq. (3.4) except that the relative sign between 
the two terms is positive. We can decompose m-^ and m 2 ^ into a pole and an non-pole 
piece. We exhibit the pole part below. For m 2 ^ we have 

m ? = -2x/D + m\ (NP) (3.5) 

where x = Y t A 2 R F/E. Similarly for m 2 ^ we have 

m| = -x/D + m| (NP) (3.6) 



From the above it is seen that although is numerically smaller than m~ 2 the pole 
part is larger in magnitude for than for m-^. Consequently the variation of is 
always larger than the variation of m~ as we approach the quasi fixed point. Further the 
proximity of the fixed point makes the transition of to the tachyonic mode rapid and 
one expects large regions of parameter space to be eliminated when one is near the fixed 
point. 

We note that the first two generations of squarks and the three generations of sleptons 
have in general a negligible correction due to the fact that the -f^-Higgs coupling to their 
corresponding quark and lepton masses is negligible. Thus for these states one does not 
expect any rapid variations in the vicinity of the Landau pole. However, rapid variations 
occur for some of the chargino and the neutralino spectrum as we approach the Landau 
pole. This can be seen most easily in the scaling limit where the chargino spectrum is 
given by 

m Wi ~ m 2 , m W2 ~ n (3.7) 

where mi = (a2/«G)^i/2- From the above we infer that will show no variation 

while rriw 2 will show a significant variation as we approach the Landau pole. Similarly the 
neutralino masses in the scaling limit are given by 

m Zl - m wJ 2 i m z 2 - m w 1 ( 3 - 8a ) 
m^ 3 ~ /I, niz 4 ~ n (3.86) 

Here again as we approach the fixed point we do not expect any rapid variation in mg 
and m^ 2 . However, m^ 3 and m^ 4 are expected to show a rapid variation in the vicinity of 
the Landau pole. These expectations are borne out in numerical analyses of the chargino 
and neutralino mass spectra. A similar analysis can be carried out for the Higgs. For 
the lightest Higgs the parameter \i enters only at the loop level and its dependence in 
the lightest Higgs mass is suppressed. Thus no rapid variation is expected for the lightest 
Higgs in the vicinity of the Landau pole. The three remaining Higgs in the scaling limit are 
all approximately degenerate with a common mass m 2 A ~ 2m 2 /sin2(3.0ne can determine 
m| from the radiative breaking equation 

sin2/3 = 2m\j{2y? + m 2 Hl + m% 2 ) (3.9a) 
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Using Eqs(2.15) and (2.16) one finds 

<=-Jm Y ^ Al k +m ^ (NP) (3 - 96) 

We see that the three heavy Higgs possess a pole piece and thus show a rapid variation 
as we approach the Landau pole. These expectations are again borne out by numerical 
computations of the mass spectra. m\ is determined by the renormalization group to be 

ml = -Boy 2 + rfj m 1 /2 + sA yo (3.10) 

Here y? = yl q and r, s, q are defined in [16]. One finds then that So becomes singular as 
we approach the Landau pole. 

Next we discuss the effect on electroweak physics due to small changes in the input 
parameters when one is in the vicinity of the Landau pole. The input parameters for the 
case under discussion are pi = (mo, m-i/2, A t , tan/3, aa, Mq, mt), i = 1, ...,7. Let us 
denote by Q a the output quantities which exhibit a pole structure near the fixed point: 

Q a = C a /D + Q a (NP) (3.11) 

We denote by Q ayi the derivative of Q a with respect to p^. Then 

Q a ^ = C a4 /D - C a D^/Dl + Q a (NP) ti (3.12) 

It is now seen that the largest variations arise for those parameters for which D ,i are non- 
vanishing. There are three such parameters: mt, olg-, and tan/3. We shall call these order 
two parameters since their variations involve a double pole. The effect on electroweak 
physics of the variations of these parameters is important since there are currently sig- 
nificant experimental errors in the determinations of, for example, mt, and «g- These 
experimental errors can be vastly exaggerated if one is in the vicinity of the Landau pole. 
We exhibit this maginification of errors in Figl. Here the m~ mass is plotted as a function 
of mt for various assumed values of tan (3. In the example shown (for the case A t = 0.5) 
one finds that the m| mass begins to show a very rapid variations as we approach the 
Landau pole. Thus in the vicinity of the Landau pole changes in m t of the order of a 
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few GeV can lead to changes in the mass of several hundred GeV. Similarly a small 
fractional change in tan/3 can generate a similar large change in m~ mass. A similar phe- 

Z i 

nomenon occurs for variations in olq. This means that precision analyses of mass spectra 
will require a very high degree of accuracy in the parameters mt, olq if one happens to lie 
in the vicinity of the Landau pole. 

Fig 1 also shows that the stop 1 turns tachyonic as we approach the Landau pole. 
It can be seen from eq(3.5) that the residue of the pole term is negative and thus as m t 
approaches m{ the pole term cancels the non-pole term in eq(3.5) turning t\ tachyonic. We 
wish to point out the important role that the trilinear coupling A t and its sign play in 
turning t\ tachyonic. Actually the parameter which controls the approach to tachyonic 
limit is A R which is defined by eq(2.11). In the vicinity of the pole x ~ A 2 R /6. To get an 
idea of the sizes of various entries in eq(2.11a) we find that for the same parameters as 
used in eq(2.6) one has 

A R ~ A t - 0.613mg (3.13) 

One finds then that for A t positive there is a cancellation between the two terms on the 
right hand side of eq(3.13) which reduces the size of Ar and slows down the approach to 
the tachyonic limit. Thus for positive A t one will have to get relatively close to the Landau 
pole for ii to turn tachyonic. For the A t negative case there is a reinforcement between the 
two terms on the right hand side of eq(3.13). Thus in this case approach to the tachyonic 
limit is much faster. An illustration of these results is given in Fig 2. Here one sees clearly 
the strong influence that the Landau pole has in turning the stop 1 tachyonic. A more 
detailed investigation of the tachyonic condition is given in Sec 4, where the tachyonic limit 
is used to delineate the allowed parameter space of MSGM. 

The coefficient of the double pole term exhibits certain scaling laws. Thus neglecting 
the single pole and the non-pole terms one has 

Qa^/Q^j = D 0ti /Do,j (3.14) 

From the above we see that the ratio of variations is independent of the index a and hence 
also independent of all the parameters that do not enter in D . If we label the order two 
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variables m t , olq-, and sinf3 by i = 1, 2, 3, then 



Q a ,i/Qa,3 = -sinf3/m t (3.15a) 

Qa,2/Q a ,i = -(l/2)m t d(log(F/£))/d(a G ) (3.156) 

We emphasize, however, that there are important non pole corrections which one must 
include in any realistic evaluation of the variations. 

We wish to investigate next the issue of how close one may approach the Landau 
pole. First from the stability of the the potential at the GUT scale one requires that 

A 2 < 3(3ml + /4) (3.16) 

Now 

^ = fi&l + /3 2 tf /b2 (l + /?it) 3/bl (A)) 1/2 (3.17) 

We have already seen that from radiative breaking of the electro-weak symmetry one 
deduces \i 2 to have the form fi 2 = a\ + c^/A)- Using these results we can deduce the D 
dependence of Hq. One finds 

$ = a/D l /2 + b/D 3 /2 (3.18) 
Using eq. (3.18) the condition of stability of the potential at the GUT scale takes the form 

ZmlDl + aDl /2 + bD^ 2 - A 2 R /3 > (3.19) 

1 /2 

If one is close enough to the pole that powers of D higher than D ' can be neglected 
then one obtains 

Dl /2 > A 2 R /(3b) (3.20) 

The above equation gives the point of closest approach to the Landau pole. However, 
as one gets very close to the fixed point other phenomenon must be taken into account. 
These include higher order terms in the top quark Yukawa coupling equations, loop effects 
in the radiative breaking equations etc. Here we point out an interesting phenomenon 
associated with inclusion of loop effects in radiative breaking and its implication in the 
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Landau pole analysis. With inclusion of one-loop effects to the effective potential, the 
condition of eq. (2.12) is modified as follows: 

, m 2 M +E 1 -(m 2 M + E 2 )tan 2 /3 1 

I? = — ^ V^ 2 - -Ml 3.21 

P tan 2 /3-l 2 z K J 

where E 1 ' 2 are the one loop corrections. The largest contributions to E 1 ' 2 arise from the 
stops and are given in the Appendix. Using results of the Appendix in eq. (3.21) we 
find that the satisfaction of radiative electro-weak symmetry breaking now shows that the 
solution of fj 2 with the above corrections is more complicated. It has the form 

H 2 = ai + b 1 /e + cilog(die)/e (3.22) 

where e = D . From eq. (3.22), we see that inclusion of the one loop correction in the 
effective potential brings in the new feature of the log(e)/e term while the tree contribution 
is only 0(l/e). This means that as we approach the Landau pole the contribution of the 
one loop effective potential eventually becomes larger than the tree contribution due to 
the log(e) divergence. This would signal the breakdown of the perturbation theory. Thus 
perturbation theory puts a natural cutoff on how close one can approach the Landau pole. 
The point of closest approach is defined by the condition that the loop contribution in 
eq. (3.22) be smaller than the tree contribution. 

4. The Parameter Space of the Minimal Supergravity Model 

In this section we shall discuss the implications of the condition that the SUSY spec- 
trum be free of tachyons in the analysis of the radiative electroweak symmetry breaking. 
As it turns out, of the 32 SUSY particles, only the light stop may become tachyonic in 
some region of the parameters. From eq. (3.5) it is easy to understand why the light stop 
can become tachyonic. Since x > 0, the first term of Eq. (3.5) will always dominate the 
second for D small enough. Thus as one approaches the Landau pole, m~ will always 
turn tachyonic, provided eq. (3.20) is not violated. One can then ask when the light stop 
becomes tachyonic. The answer is very sensitive to the values one chooses for the param- 
eters, since c, e, /, /i, /2, fz, h, and k have very complicated dependence on Mg, ag, 
m t (via Yf), and tan/3. To determine if the light stop is tachyonic, we need to solve the 
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equation = 0. Fortunately, this equation can be solved analytically for A t = 0. We 
consider two points in the parameter space for A t = 0, tan/3 = 1.732, and m t = 150 GeV. 

Case l):We choose the remaining parameters to be Mq = 10 GeV, aa = 1/25.7 
(we may replace M§ cos 2 9w by while computing Y t at the electroweak scale), and 
«2 = 0.03322. Scaling all the masses by 100 GeV, the light stop has zero mass on the m 1( / 2 
and mo plane along the following line: 

19.66 - 15.11m 2 /2 - A.hAm\ /2 + 3.77m 2 , + 0.74m 2 /2 m^ + 0.39m£ = 0, (4.1) 

Case 2): We choose remaining parameters to be Mq = 10 16187 GeV, ao = 1/24.11, 
and «2 = 0.03358. Using the same scaling as in case 1, the light stop has zero mass on the 
mi/2 and mo plane along the following line: 

19.04 + 28.44m 2 /2 + 18Almj /2 + 6.16m 2 , + 6A7m 2 1/2 m 2 + 0.55m^ = 0, (4.2) 

There is a significant difference between Case 1 and Case 2. For Case 1, coefficients 
in eq (4.1) appear with both positive and negative signs, while in Case 2, coefficients in 
eq (4.2) all have positive signs. Thus, the condition that the light stop be tachyonic may 
be satisfied in Case 1 for points below the line given by equation (4.1), while this same 
condition cannot be satisfied at any point on the mo-mi/ 2 plane for Case 2, even though 
the two cases differ by only a small change in ac and other parameters. Cases 1 and 
2 discussed above were for illustrative purposes only. We discuss next a more systematic 
study of the condition that the light stop be nontachyonic. Figures 3 and 4 show the results 
for mt = 170 GeV (corresponding to a pole mass of about 175 GeV), and tan/3 = 5.0. 
Fig. 3 shows the results for A t /mo = —0.25, —0.3, —0.4, —0.5, and —0.8 on the m^-mo 
plane, the light stop is tachyonic above and right to each curve for a fixed A t . Fig. 4 
demonstrates the results for A t /m = 0.5, 0.8, 1.0, 2.0, and 5.0 on the mg-m plane, the 
light stop is tachyonic above and left to each curve for a fixed A t . From these figures, one 
observes that the condition that the light stop be nontachyonic excludes more parameter 
space when \x and A t have the same sign than when they have opposite signs. Further large 
values of A t of either sign, specifically values of A t < — 0.8m and A t > 5.0m , are excluded 
by this condition. Thus the allowed region of A t lies within — 0.8mo < A t < 5.0mo for 
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tanf3=5 and m t =170 GeV. This is a very strong constraint and it is interesting that the 
light stop being tachyonic alone can impose such a stringent constraint on the parameter 
space of MSGM. 



5. Conclusions 

We have investigated in this paper the phenomena associated with effects of the 
Landau pole in the top quark Yukawa coupling and implications of the constraint that the 
SUSY spectrum be free of tachyons. Regarding the first we find that in the renormalization 
group analyses some of the mass spectra (i.e., masses of the stops, the heavy chargino, 
heavy neutralinos, and the heavy Higgs) are very sensitive to the input errors in mt, 
etc, and tan/3 when the top mass is close to its fixed point value. In this region of the 
parameter space precision physics requires much greater accuracy of the input data. We 
also investigated the issue of how close one may approach the fixed point. Regarding the 
second phenomenon the condition that the SUSY spectrum be tachyon free turns out to 
be a rather stringent one as it eliminates a large part of the parameter space. For example, 
the trilinear scalar coupling, i.e., A t is confined within the range of — 0.8mo to 5.0mo, for 
the parameters of Fig. 3 and 4. The graphs indicate the strong correlation between At 
and //, i.e., the light stop is more likely to become tachyonic when A t and \i have the same 
sign. 
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Appendix 

The contribution to E 1 from stops £i 5 2 is given by 
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Similarly the contribution to E 2 from stops is 

«2 + «Y 9 r~ 1 777.? 1 r l/9 9wl 4 ., 

£? =3-^- — -m, 2 — + 9 9 =F t — 9 ^{r m - m 2 )(-- + - sin 2 W ) 

4l . a 8tt 4 M 2 sin 2 /3 m 2 -m? l 2 v *« A 2 3 ; 

777^ 777^- 

+ ^2 2, m o A f (/x cot /3 + m A t )}] log(-^). (A2) 
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Table Captions 

Table 1 The fit to the electroweak data for the gauge coupling constants given in equa- 
tion (2.9) in one and two loop approximations, for Mq = 10 16187 GeV, and various values 
for da- 

Table 2 The effects of the variation of cxgut on the ratio F/E and hence on the top 
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Landau pole in one and two loop R.G. analysis. One finds that the variations of oigut are 
much more significant than the two loops effects. 



Figure Captions 



Fig. 1 Exhibition of m~ vs.mt for the case when A t = 0.5mo, mo = 600 GeV, mg = 300 
GeV and tan/? = 1.2, 1.4, 2.0, 3.0, 5.0, 7.0, 9.0 in increasing value as we go left to right. We 
see that there is a huge magnification of error in mj for a corresponding small error in m t 

^1 

when the top mass is close to its fixed point value. 



Fig. 2 Exhibition of the effect of the Landau pole on the tachyonic limit for the case 
m =600 GeV,m g i u i no =300GeV and \i < 0. The top curve gives the position of the 
Landau pole as a function of tan(3. The middle curve gives, for the case when A t = 0.5mo, 
the value of mt as a function of tan(3 so that for all mt values above this value m~ is always 
tachyonic. The bottom curve is the same as the middle curve except that A t = — 0.5mo- 



Fig. 3 The solutions for m| = on the m d -m Q plane for a G = 24.11, M G = 10 16187 
GeV, and mt = 170 GeV and tan/? = 5.0. This figure shows the branch for At > 0. The 
light stop is tachyonic above and left to each curve for a fixed A t . All masses are in GeV. 



Fig. 4 The solutions for m~ = on the m„-mo plane for a G = 24.11, Mq = io 16187 

ti y 

GeV, and m t = 170 GeV and tan/? = 5.0. This figure shows the branch for A t < 0. The 
light stop is tachyonic above and right to each curve for a fixed A t . All masses are in GeV. 
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Table 1 





One Loop 


Two Loop 


(Xq 1 = 25.7 


0.0166359 
0.0323481 
0.0994 


0.0165156 
0.031735 


a^ 1 = 24.5 

OL\ 

a 3 


0.0169748 
0.0336544 
0.1129 


0.0168391 
0.0329366 


a^ 1 = 24.1 

a 2 


0.0170909 
0.0341137 
0.1182 


0.0169493 
0.0333552 


a" 1 = 24.0 

«2 
«3 


0.0171201 
0.0342304 
0.1196 


0.016977 
0.0334613 



Table 2 





One Loop 


Two Loop 


(Xq 1 = 24.11 






E 


12.985 


13.290 


F 


266.48 


275.97 


F/E 


20.523 


20.766 


cxq 1 = 24.5 






E 


11.938 


12.536 


F 


252.61 


266.78 


F/E 


21.160 


21.288 


a^ 1 = 25.0 






E 


11.101 


12.536 


F 


242.46 


256.15 


F/E 


21.842 


21.919 
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